\documentclass{article}
\usepackage[utf8]{inputenc}
\usepackage[UKenglish]{babel}
\usepackage{%
xcolor,%
framed,%
xspace,%
mdframed,%
url%
}

\definecolor{light-gray}{gray}{0.90}
\definecolor{shadecolor}{rgb}{1,0.9,0.7}

\title{Analysing termite favouritism in annular arenas\footnote{This is part of the paper: Paiva \textit{et al.}(in review) Scale-free movement patterns in termites emerge from social interactions and preferential attachments. \textit{submitted to PNAS}. }}



\usepackage{Sweave}
\begin{document}

\input{graphs/fig-concordance}
\maketitle
\begin{description}
\item[Compiled data and ran the analysis:]Og DeSouza \url{<og.souza@ufv.br>}
\item[Tracking algorithm:]\ 
\begin{description}
\item[Available at:]\url{https://github.com/dmrib/trackingtermites} 
\item[Wrote the final code:] Danilo Ribeiro
\item[Wrote the first working code:] Sidiney Alves
\item[Supervised coding:] Og DeSouza
\end{description}
\end{description}
\section{Overall rationale}
%\url{https://github.com/dmrib/trackingtermites}

\begin{itemize}
\item Here we inspect whether termites walking in annular arenas would preferentially contact some of their nestmates over others along a \textit{c.a.} 30 min period.
\item our tracking algorithm (available at GitHub \url{https://github.com/dmrib/trackingtermites}) was able to (i) identify each one of the termite individuals in the arena, (ii) track their positions along the whole footage, and (iii) tally the number of time-steps each individual termite acted as a target of a one-to-one interaction. That is to say, when a given termite walked towards another and met it, this second termite was recorded as a ``target". We then recorded the duration of this encounter in terms of time-steps. Then we analysed the proportion of time-steps (relative to the duration of the whole footage) \textit{each} termite acted as a target. In the absence of favouritism, the proportion of time-steps each member of the group acted as a target would not differ among group members (that is our null hypothesis). That is, for a given footage duration, target time-steps would be even distributed among all members of the group. On the other hand, if *some* of the termites in the group would be "preferred targets" over others, they would spend more time-steps acting as targets than other members of the group. In doing so, for a given footage duration, target time-steps would be unevenly distributed among all members of the group. 

\item To do so, we filmed and tracked \underline{each individual} termite in the confined group, tallying the number of time-steps this focal termite spent contacting a given target termite along the whole footage. 
\item This has produced a dataset such as:
\begin{verbatim}
Tallyied time-steps for a group of 3 individuals.
--------------------------------------------------
  trajectory focalTermite targetTermite steps
1  traj00006    termite01     termite02   640
2  traj00006    termite01     termite03   676
3  traj00006    termite02     termite01   635
4  traj00006    termite02     termite03   500
5  traj00006    termite03     termite01   702
6  traj00006    termite03     termite02   502
--------------------------------------------------
\end{verbatim}
\item then, for each individual we averaged number of time-steps it acted as a \underline{target} of an interaction along the whole footage, producing a result such as:
\begin{verbatim}
Averaged time-steps a termite spent as a target.
------------------------------------------------
  targetTermite steps
1     termite01 668.5
2     termite02 571.0
3     termite03 588.0
------------------------------------------------

\end{verbatim}
\item then we ran a Chi-square test to verify whether the proportion of time-steps spent as a target would vary among termite individuals. That is to say, we inspected whether these proportions would depart from an uniform distribution and, if so, to point to the existence of `favourite targets'.
\item in case the Chi-square test did point to significant differences, we counted the number of `favourites' within that given group. These `favourites' were defined as as being the termites whose average number of time-steps spent as a target was larger than the global average of targeting time-steps in the whole footage.
\item The number of preferred target-termites in each arena (y-var) was then regressed against the number of termites therein confined (x-var), in order to inspect how preferential attachments would correlate with group size. Analysis consisted in contrasting the model thereby obtained with a model with zero intercept and slope $=$ 1, that is, a model in which favouritism was absent. Modelling was performed in R, under Generalised Linear Modelling and normal error distribution, followed by residual analysis. Contrasts were made using Akaike Information Criterion (AIC).

\end{itemize}
% We inspected the frequency of interindividual contacts between termites confined together in an annular arena, to verify whether they would favour a given conspecific over the other members of the group.
% 
% To to so, we filmed the termite group over, in average, 34m12s (min $=$ 20m52s; max $=$ 53m12s). We then tracked each individual in the group, tallying the number of time-steps it spent contacting each one of the others. A Python algorithm based on OpenCV was customised to track the termite identities and respective time-steps a given focal termite spent contacting a given target termite. This was done for ALL termites in the group.
% 
% The analysis proceeded inspecting whether the total amount of time-steps a given \textit{target termite} was 
% 
% 
% proportion of time-steps a given focal-target pair of termites spent together would divert from an equal share




\section{Loading and preparing data}
\begin{Schunk}
\begin{Sinput}
> setwd("./")
> dt<-read.table("../dat/newdata.dat",sep=",") # termite1 changed to termite01, etc
> names(dt)[1] <- "trajectory"
> names(dt)[2]<- "focalTermite"
> names(dt)[3]<- "targetTermite"
> names(dt)[4]<- "steps"
> summary(dt)
\end{Sinput}
\begin{Soutput}
     trajectory      focalTermite    targetTermite     steps       
 traj00013:240    termite01: 67    termite01: 67   Min.   :   1.0  
 traj00130:240    termite02: 67    termite02: 67   1st Qu.: 309.0  
 traj00007: 56    termite03: 67    termite03: 67   Median : 626.0  
 traj00100: 56    termite04: 63    termite04: 63   Mean   : 858.4  
 traj00145: 56    termite05: 51    termite05: 51   3rd Qu.:1032.0  
 traj00012: 12    termite06: 51    termite06: 51   Max.   :8963.0  
 (Other)  : 48   (Other)   :342   (Other)   :342                   
\end{Soutput}
\end{Schunk}
\newpage
\section{Preferences per trajectory}
Preparing functions:

\begin{Schunk}
\begin{Sinput}
> interact <- function(traj){ # interactions with target termite
+   ss <- subset(dt,trajectory==traj)
+   ss <- droplevels(ss)
+   ##
+   ## the total number of steps here considers readings
+   ## from termite01 > termite02 as well as termite02 > termite01.
+   ## We should make it as the mean of both cases.
+   totSteps <- sum(ss$steps)
+   interactions <- aggregate(ss$steps,list(ss$targetTermite),mean)
+   names(interactions)[1] <- "targetTermite"
+   names(interactions)[2] <- "steps"
+   return(interactions)
+ }
> pl<-function(interactions)
+   {pval <- round(chisq.test(interactions[2])$p.value,4)
+   expected <- chisq.test(interactions[2])$expected
+   barplot(interactions$steps,
+         names.arg = interactions$targetTermite,
+         las=2,
+         main=paste(traj, " p = ",pval ),
+         ylab="time-steps",
+         #xlab="termite ID",
+         ylim=c(0,1.1*max(interactions[2]))
+         )
+ abline(h=expected[1],col="red")
+ }
> summtraj <- function(resultChi){
+   expected <- resultChi$expected
+   observed <- resultChi$observed
+   favourites <- sum(observed>expected)
+   total <- length(resultChi$observed)
+   totSteps <- sum(resultChi$observed)
+   chi <- resultChi$statistic
+   pval <- resultChi$p.value
+   return(list(total = total, 
+               favourites = favourites, 
+               totSteps = totSteps,
+               chisquare = chi[[1]], 
+               pvalue = pval))
+ }
> 
\end{Sinput}
\end{Schunk}
\newpage
\section{Trajectory 00006}
\begin{Schunk}
\begin{Sinput}
> traj <- "traj00006"
> ii <- interact(traj)
> ii
\end{Sinput}
\begin{Soutput}
  targetTermite steps
1     termite01 668.5
2     termite02 571.0
3     termite03 588.0
\end{Soutput}
\begin{Sinput}
> resultChi <- chisq.test(ii[2])
> resultTraj00006 <-summtraj(resultChi)
> resultTraj00006
\end{Sinput}
\begin{Soutput}
$total
[1] 3

$favourites
[1] 1

$totSteps
[1] 1827.5

$chisquare
[1] 8.905882

$pvalue
[1] 0.01164427
\end{Soutput}
\begin{Sinput}
> pl(ii)
> svg(filename = "graphs/fig-traj00006.svg")
> pl(ii)
> dev.off()
\end{Sinput}
\begin{Soutput}
pdf 
  2 
\end{Soutput}
\end{Schunk}
\includegraphics{graphs/fig-traj00006}

\section{Trajectory 00007}
\begin{Schunk}
\begin{Sinput}
> traj <- "traj00007"
> ii <- interact(traj)
> resultChi <- chisq.test(ii[2])
> resultTraj00007 <-summtraj(resultChi)
> resultTraj00007
\end{Sinput}
\begin{Soutput}
$total
[1] 8

$favourites
[1] 3

$totSteps
[1] 2482.571

$chisquare
[1] 131.5284

$pvalue
[1] 3.012616e-25
\end{Soutput}
\begin{Sinput}
> pl(ii)
> svg(filename = "graphs/fig-traj00007.svg")
> pl(ii)
> dev.off()
\end{Sinput}
\begin{Soutput}
pdf 
  2 
\end{Soutput}
\end{Schunk}
\includegraphics{graphs/fig-traj00007}


\section{Trajectory 00012}
\begin{Schunk}
\begin{Sinput}
> traj <- "traj00012"
> ii <- interact(traj)
> resultChi <- chisq.test(ii[2])
> resultTraj00012 <-summtraj(resultChi)
> resultTraj00012
\end{Sinput}
\begin{Soutput}
$total
[1] 4

$favourites
[1] 2

$totSteps
[1] 3315.667

$chisquare
[1] 36.35384

$pvalue
[1] 6.303287e-08
\end{Soutput}
\begin{Sinput}
> pl(ii)
> svg(filename = "graphs/fig-traj00012.svg")
> pl(ii)
> dev.off()
\end{Sinput}
\begin{Soutput}
pdf 
  2 
\end{Soutput}
\end{Schunk}
\includegraphics{graphs/fig-traj00012}

\section{Trajectory 00013}
\begin{Schunk}
\begin{Sinput}
> traj <- "traj00013"
> ii <- interact(traj)
> resultChi <- chisq.test(ii[2])
> resultTraj00013 <-summtraj(resultChi)
> resultTraj00013
\end{Sinput}
\begin{Soutput}
$total
[1] 16

$favourites
[1] 9

$totSteps
[1] 4806

$chisquare
[1] 46.02939

$pvalue
[1] 5.260081e-05
\end{Soutput}
\begin{Sinput}
> pl(ii)
> svg(filename = "graphs/fig-traj00013.svg")
> pl(ii)
> dev.off()
\end{Sinput}
\begin{Soutput}
pdf 
  2 
\end{Soutput}
\end{Schunk}
\includegraphics{graphs/fig-traj00013}

\section{Trajectory 00095}
\begin{Schunk}
\begin{Sinput}
> traj <- "traj00095"
> ii <- interact(traj)
> resultChi <- chisq.test(ii[2])
> resultTraj00095 <-summtraj(resultChi)
> resultTraj00095
\end{Sinput}
\begin{Soutput}
$total
[1] 4

$favourites
[1] 2

$totSteps
[1] 5805

$chisquare
[1] 72.91237

$pvalue
[1] 1.014982e-15
\end{Soutput}
\begin{Sinput}
> pl(ii)
> svg(filename = "graphs/fig-traj00095.svg")
> pl(ii)
> dev.off()
\end{Sinput}
\begin{Soutput}
pdf 
  2 
\end{Soutput}
\end{Schunk}
\includegraphics{graphs/fig-traj00095}

\section{Trajectory 00096}
\begin{Schunk}
\begin{Sinput}
> traj <- "traj00096"
> ii <- interact(traj)
> resultChi <- chisq.test(ii[2])
> resultTraj00096 <-summtraj(resultChi)
> resultTraj00096
\end{Sinput}
\begin{Soutput}
$total
[1] 4

$favourites
[1] 1

$totSteps
[1] 3808

$chisquare
[1] 103.4799

$pvalue
[1] 2.774122e-22
\end{Soutput}
\begin{Sinput}
> pl(ii)
> svg(filename = "graphs/fig-traj00096.svg")
> pl(ii)
> dev.off()
\end{Sinput}
\begin{Soutput}
pdf 
  2 
\end{Soutput}
\end{Schunk}
\includegraphics{graphs/fig-traj00096}

\section{Trajectory 00100}
\begin{Schunk}
\begin{Sinput}
> traj <- "traj00100"
> ii <- interact(traj)
> resultChi <- chisq.test(ii[2])
> resultTraj00100 <-summtraj(resultChi)
> resultTraj00100
\end{Sinput}
\begin{Soutput}
$total
[1] 8

$favourites
[1] 4

$totSteps
[1] 22134.86

$chisquare
[1] 1149.136

$pvalue
[1] 7.031052e-244
\end{Soutput}
\begin{Sinput}
> pl(ii)
> svg(filename = "graphs/fig-traj00100.svg")
> pl(ii)
> dev.off()
\end{Sinput}
\begin{Soutput}
pdf 
  2 
\end{Soutput}
\end{Schunk}
\includegraphics{graphs/fig-traj00100}

\section{Trajectory 00130}
\begin{Schunk}
\begin{Sinput}
> traj <- "traj00130"
> ii <- interact(traj)
> resultChi <- chisq.test(ii[2])
> resultTraj00130 <-summtraj(resultChi)
> resultTraj00130
\end{Sinput}
\begin{Soutput}
$total
[1] 16

$favourites
[1] 7

$totSteps
[1] 15451.27

$chisquare
[1] 93.37328

$pvalue
[1] 2.319141e-13
\end{Soutput}
\begin{Sinput}
> pl(ii)
> svg(filename = "graphs/fig-traj00130.svg")
> pl(ii)
> dev.off()
\end{Sinput}
\begin{Soutput}
pdf 
  2 
\end{Soutput}
\end{Schunk}
\includegraphics{graphs/fig-traj00130}

\section{Trajectory 00143}
\begin{Schunk}
\begin{Sinput}
> traj <- "traj00143"
> ii <- interact(traj)
> resultChi <- chisq.test(ii[2])
> resultTraj00143 <-summtraj(resultChi)
> resultTraj00143
\end{Sinput}
\begin{Soutput}
$total
[1] 3

$favourites
[1] 2

$totSteps
[1] 4616

$chisquare
[1] 118.3383

$pvalue
[1] 2.009875e-26
\end{Soutput}
\begin{Sinput}
> pl(ii)
> svg(filename = "graphs/fig-traj00143.svg")
> pl(ii)
> dev.off()
\end{Sinput}
\begin{Soutput}
pdf 
  2 
\end{Soutput}
\end{Schunk}
\includegraphics{graphs/fig-traj00143}

\section{Trajectory 00144}
\begin{Schunk}
\begin{Sinput}
> traj <- "traj00144"
> ii <- interact(traj)
> resultChi <- chisq.test(ii[2])
> resultTraj00144 <-summtraj(resultChi)
> resultTraj00144
\end{Sinput}
\begin{Soutput}
$total
[1] 4

$favourites
[1] 3

$totSteps
[1] 17279.33

$chisquare
[1] 3145.936

$pvalue
[1] 0
\end{Soutput}
\begin{Sinput}
> pl(ii)
> svg(filename = "graphs/fig-traj00144.svg")
> pl(ii)
> dev.off()
\end{Sinput}
\begin{Soutput}
pdf 
  2 
\end{Soutput}
\end{Schunk}
\includegraphics{graphs/fig-traj00144}

\section{Trajectory 00145}
\begin{Schunk}
\begin{Sinput}
> traj <- "traj00145"
> ii <- interact(traj)
> resultChi <- chisq.test(ii[2])
> resultTraj00145 <-summtraj(resultChi)
> resultTraj00145
\end{Sinput}
\begin{Soutput}
$total
[1] 8

$favourites
[1] 4

$totSteps
[1] 4012.714

$chisquare
[1] 47.49812

$pvalue
[1] 4.462514e-08
\end{Soutput}
\begin{Sinput}
> pl(ii)
> svg(filename = "graphs/fig-traj00145.svg")
> pl(ii)
> dev.off()
\end{Sinput}
\begin{Soutput}
pdf 
  2 
\end{Soutput}
\end{Schunk}
\includegraphics{graphs/fig-traj00145}

\section{Big summary}
\begin{Schunk}
\begin{Sinput}
> traj <- sort(unique(dt$trajectory))
> groupSize <- c(resultTraj00006$total,
+                resultTraj00007$total,
+                resultTraj00012$total,
+                resultTraj00013$total,
+                resultTraj00095$total,
+                resultTraj00096$total,
+                resultTraj00100$total,
+                resultTraj00130$total,
+                resultTraj00143$total,
+                resultTraj00144$total,
+                resultTraj00145$total
+                )
> favourites <- c(resultTraj00006$favourites,
+                 resultTraj00007$favourites,
+                 resultTraj00012$favourites,
+                 resultTraj00013$favourites,
+                 resultTraj00095$favourites,
+                 resultTraj00096$favourites,
+                 resultTraj00100$favourites,
+                 resultTraj00130$favourites,
+                 resultTraj00143$favourites,
+                 resultTraj00144$favourites,
+                 resultTraj00145$favourites
+                 )
> steps <- c(resultTraj00006$totSteps,
+                 resultTraj00007$totSteps,
+                 resultTraj00012$totSteps,
+                 resultTraj00013$totSteps,
+                 resultTraj00095$totSteps,
+                 resultTraj00096$totSteps,
+                 resultTraj00100$totSteps,
+                 resultTraj00130$totSteps,
+                 resultTraj00143$totSteps,
+                 resultTraj00144$totSteps,
+                 resultTraj00145$totSteps
+                 )
> chisq <- c(resultTraj00006$chisquare,
+                 resultTraj00007$chisquare,
+                 resultTraj00012$chisquare,
+                 resultTraj00013$chisquare,
+                 resultTraj00095$chisquare,
+                 resultTraj00096$chisquare,
+                 resultTraj00100$chisquare,
+                 resultTraj00130$chisquare,
+                 resultTraj00143$chisquare,
+                 resultTraj00144$chisquare,
+                 resultTraj00145$chisquare
+                 )
> pval <- c(resultTraj00006$pvalue,
+                 resultTraj00007$pvalue,
+                 resultTraj00012$pvalue,
+                 resultTraj00013$pvalue,
+                 resultTraj00095$pvalue,
+                 resultTraj00096$pvalue,
+                 resultTraj00100$pvalue,
+                 resultTraj00130$pvalue,
+                 resultTraj00143$pvalue,
+                 resultTraj00144$pvalue,
+                 resultTraj00145$pvalue
+                 )
> summdf <- data.frame(traj,groupSize,fav=favourites, steps=round(steps,0), chisq, pval)          
> nestID <- read.table("nestIDtraj.csv", h=T)
> summTab <- merge(nestID,summdf,by="traj")
> library(xtable)
> print(xtable(summTab,
+              display = c("s","s","s","s","s","s","g","g")),
+       math.style.exponents=TRUE
+       )
\end{Sinput}
% latex table generated in R 3.6.3 by xtable 1.8-4 package
% Tue Feb  9 08:28:15 2021
\begin{table}[ht]
\centering
\begin{tabular}{rllrrrrr}
  \hline
 & traj & nest & groupSize & fav & steps & chisq & pval \\ 
  \hline
1 & traj00006 & nest2017-05-24-KRR-N01 &  3 & 1 &  1828 & 8.9 & 0.012 \\ 
  2 & traj00007 & nest2017-05-24-KRR-N01 &  8 & 3 &  2483 & $1.3 \times 10^{2}$ & $3 \times 10^{-25}$ \\ 
  3 & traj00012 & nest2017-05-26-LPNS-N01 &  4 & 2 &  3316 &  36 & $6.3 \times 10^{-8}$ \\ 
  4 & traj00013 & nest2017-05-29-KRR-N01 & 16 & 9 &  4806 &  46 & $5.3 \times 10^{-5}$ \\ 
  5 & traj00095 & nest2017-05-03-KRR-N01 &  4 & 2 &  5805 &  73 & $1 \times 10^{-15}$ \\ 
  6 & traj00096 & nest2017-05-03-KRR-N01 &  4 & 1 &  3808 & $1 \times 10^{2}$ & $2.8 \times 10^{-22}$ \\ 
  7 & traj00100 & nest2017-05-03-KRR-N01 &  8 & 4 & 22135 & $1.1 \times 10^{3}$ & $7 \times 10^{-244}$ \\ 
  8 & traj00130 & nest2017-05-08-KRR-N02 & 16 & 7 & 15451 &  93 & $2.3 \times 10^{-13}$ \\ 
  9 & traj00143 & nest207-06-06-LPNS-N03 &  3 & 2 &  4616 & $1.2 \times 10^{2}$ & $2 \times 10^{-26}$ \\ 
  10 & traj00144 & nest2017-06-07-KRR-N01 &  4 & 3 & 17279 & $3.1 \times 10^{3}$ &   0 \\ 
  11 & traj00145 & nest2016-06-07-KRR-N01 &  8 & 4 &  4013 &  47 & $4.5 \times 10^{-8}$ \\ 
   \hline
\end{tabular}
\end{table}\begin{Sinput}
> write.csv(summTab, file = "summaryFavourites.csv", row.names = FALSE)
> favourites/groupSize
\end{Sinput}
 [1] 0.3333333 0.3750000 0.5000000 0.5625000 0.5000000 0.2500000 0.5000000
 [8] 0.4375000 0.6666667 0.7500000 0.5000000\begin{Sinput}
> mean(favourites/groupSize)
\end{Sinput}
[1] 0.4886364\end{Schunk}

\section{Regression}
\begin{Schunk}
\begin{Sinput}
> m1<-glm(favourites~groupSize)
> m0<-glm(favourites~offset(0+groupSize))
> AIC(m0,m1)
\end{Sinput}
\begin{Soutput}
   df      AIC
m0  2 54.68087
m1  3 29.40147
\end{Soutput}
\begin{Sinput}
> anova(m0,m1, test="F")
\end{Sinput}
\begin{Soutput}
Analysis of Deviance Table

Model 1: favourites ~ offset(0 + groupSize)
Model 2: favourites ~ groupSize
  Resid. Df Resid. Dev Df Deviance      F    Pr(>F)    
1        10     64.545                                 
2         9      5.406  1    59.14 98.466 3.816e-06 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
\end{Soutput}
\begin{Sinput}
> pfav <- anova(m0,m1, test="F")$'Pr(>F)'[2]
> library(scales)
> # digits = 3 is the default but I am setting it here to be explicit,
> # and draw attention to the fact this is different than the formatC
> # solution.
> pfav <- scientific(pfav, digits = 3)
> pfav
\end{Sinput}
\begin{Soutput}
[1] "3.82e-06"
\end{Soutput}
\begin{Sinput}
> anova(m1,test="F")
\end{Sinput}
\begin{Soutput}
Analysis of Deviance Table

Model: gaussian, link: identity

Response: favourites

Terms added sequentially (first to last)


          Df Deviance Resid. Df Resid. Dev      F    Pr(>F)    
NULL                         10     62.727                     
groupSize  1   57.322         9      5.406 95.438 4.344e-06 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
\end{Soutput}
\begin{Sinput}
> coef(m1)
\end{Sinput}
\begin{Soutput}
(Intercept)   groupSize 
-0.06323185  0.49609680 
\end{Soutput}
\begin{Sinput}
> plot(favourites~jitter(groupSize,factor=2)
+      ,lwd = 3
+      ,cex = 3
+      ,pch = 21
+      ,xlab = "number of termites in arena"
+      ,ylab = "number of favourites"
+      )
> curve(coef(m1)[1]+coef(m1)[2]*x, add=T,lwd=3) 
> curve(0+1*x, add = T, col = "red", lwd = 3, lty=2)
> eq0 <- c('y ~ x')
> eq0
\end{Sinput}
\begin{Soutput}
[1] "y ~ x"
\end{Soutput}
\begin{Sinput}
> text(5.7,7,paste("no favouritism:", eq0,sep='\n'),col="red", adj=1)
> eq1 <- paste('y ~'
+     ,round(coef(m1)[1],3)
+     ,'+'
+     ,round(coef(m1)[2],3)
+     , 'x'
+     )
> eq1
\end{Sinput}
\begin{Soutput}
[1] "y ~ -0.063 + 0.496 x"
\end{Soutput}
\begin{Sinput}
> pval <- paste('P = ', pfav)
> text(12,4, paste(
+     'favouritism:'
+     ,eq1
+     ,pval
+     ,sep='\n'
+   )
+   )
\end{Sinput}
\end{Schunk}
\includegraphics{graphs/fig-favouritismLine}

\section*{Raw data}
\begin{Schunk}
\begin{Sinput}
> dt
\end{Sinput}
\begin{Soutput}
    trajectory focalTermite targetTermite steps
1    traj00006    termite01     termite02   640
2    traj00006    termite01     termite03   676
3    traj00006    termite02     termite01   635
4    traj00006    termite02     termite03   500
5    traj00006    termite03     termite01   702
6    traj00006    termite03     termite02   502
7    traj00007    termite01     termite02   538
8    traj00007    termite01     termite03   262
9    traj00007    termite01     termite04   166
10   traj00007    termite01     termite05   306
11   traj00007    termite01     termite06   142
12   traj00007    termite01     termite07   286
13   traj00007    termite01     termite08   319
14   traj00007    termite02     termite01   535
15   traj00007    termite02     termite03   623
16   traj00007    termite02     termite04    93
17   traj00007    termite02     termite05   284
18   traj00007    termite02     termite06   452
19   traj00007    termite02     termite07   405
20   traj00007    termite02     termite08   301
21   traj00007    termite03     termite01   260
22   traj00007    termite03     termite02   630
23   traj00007    termite03     termite04   165
24   traj00007    termite03     termite05   531
25   traj00007    termite03     termite06   462
26   traj00007    termite03     termite07   538
27   traj00007    termite03     termite08   334
28   traj00007    termite04     termite01   157
29   traj00007    termite04     termite02    92
30   traj00007    termite04     termite03   166
31   traj00007    termite04     termite05    59
32   traj00007    termite04     termite06    49
33   traj00007    termite04     termite07   452
34   traj00007    termite04     termite08   152
35   traj00007    termite05     termite01   306
36   traj00007    termite05     termite02   276
37   traj00007    termite05     termite03   533
38   traj00007    termite05     termite04    60
39   traj00007    termite05     termite06   408
40   traj00007    termite05     termite07   196
41   traj00007    termite05     termite08   274
42   traj00007    termite06     termite01   150
43   traj00007    termite06     termite02   449
44   traj00007    termite06     termite03   479
45   traj00007    termite06     termite04    48
46   traj00007    termite06     termite05   428
47   traj00007    termite06     termite07   176
48   traj00007    termite06     termite08   401
49   traj00007    termite07     termite01   293
50   traj00007    termite07     termite02   399
51   traj00007    termite07     termite03   533
52   traj00007    termite07     termite04   453
53   traj00007    termite07     termite05   195
54   traj00007    termite07     termite06   168
55   traj00007    termite07     termite08   319
56   traj00007    termite08     termite01   309
57   traj00007    termite08     termite02   297
58   traj00007    termite08     termite03   333
59   traj00007    termite08     termite04   151
60   traj00007    termite08     termite05   262
61   traj00007    termite08     termite06   397
62   traj00007    termite08     termite07   326
63   traj00012    termite01     termite02   394
64   traj00012    termite01     termite03   630
65   traj00012    termite01     termite04  1084
66   traj00012    termite02     termite01   403
67   traj00012    termite02     termite03  1118
68   traj00012    termite02     termite04   808
69   traj00012    termite03     termite01   633
70   traj00012    termite03     termite02  1122
71   traj00012    termite03     termite04   910
72   traj00012    termite04     termite01  1099
73   traj00012    termite04     termite02   835
74   traj00012    termite04     termite03   911
75   traj00095    termite01     termite02  1765
76   traj00095    termite01     termite03  1500
77   traj00095    termite01     termite04  1024
78   traj00095    termite02     termite01  1744
79   traj00095    termite02     termite03  1851
80   traj00095    termite02     termite04  1238
81   traj00095    termite03     termite01  1505
82   traj00095    termite03     termite02  1894
83   traj00095    termite03     termite04  1357
84   traj00095    termite04     termite01   978
85   traj00095    termite04     termite02  1238
86   traj00095    termite04     termite03  1321
87   traj00096    termite01     termite02   537
88   traj00096    termite01     termite03  1040
89   traj00096    termite01     termite04  1084
90   traj00096    termite02     termite01   536
91   traj00096    termite02     termite03   489
92   traj00096    termite02     termite04  1540
93   traj00096    termite03     termite01  1028
94   traj00096    termite03     termite02   497
95   traj00096    termite03     termite04  1046
96   traj00096    termite04     termite01  1064
97   traj00096    termite04     termite02  1516
98   traj00096    termite04     termite03  1047
99   traj00100    termite01     termite02  2190
100  traj00100    termite01     termite03  3335
101  traj00100    termite01     termite04  3485
102  traj00100    termite01     termite05  3294
103  traj00100    termite01     termite06  2317
104  traj00100    termite01     termite07  5373
105  traj00100    termite01     termite08  1822
106  traj00100    termite02     termite01  2206
107  traj00100    termite02     termite03  2266
108  traj00100    termite02     termite04  1910
109  traj00100    termite02     termite05  3978
110  traj00100    termite02     termite06  2914
111  traj00100    termite02     termite07  2262
112  traj00100    termite02     termite08  1931
113  traj00100    termite03     termite01  3338
114  traj00100    termite03     termite02  2265
115  traj00100    termite03     termite04  4718
116  traj00100    termite03     termite05  6029
117  traj00100    termite03     termite06  1220
118  traj00100    termite03     termite07  2186
119  traj00100    termite03     termite08  1392
120  traj00100    termite04     termite01  3497
121  traj00100    termite04     termite02  1891
122  traj00100    termite04     termite03  4725
123  traj00100    termite04     termite05  3877
124  traj00100    termite04     termite06  1470
125  traj00100    termite04     termite07  1684
126  traj00100    termite04     termite08  1770
127  traj00100    termite05     termite01  3298
128  traj00100    termite05     termite02  3956
129  traj00100    termite05     termite03  6030
130  traj00100    termite05     termite04  3860
131  traj00100    termite05     termite06  4038
132  traj00100    termite05     termite07  4877
133  traj00100    termite05     termite08  1899
134  traj00100    termite06     termite01  2317
135  traj00100    termite06     termite02  2909
136  traj00100    termite06     termite03  1211
137  traj00100    termite06     termite04  1470
138  traj00100    termite06     termite05  4034
139  traj00100    termite06     termite07  1614
140  traj00100    termite06     termite08  1348
141  traj00100    termite07     termite01  5350
142  traj00100    termite07     termite02  2258
143  traj00100    termite07     termite03  2180
144  traj00100    termite07     termite04  1686
145  traj00100    termite07     termite05  4866
146  traj00100    termite07     termite06  1614
147  traj00100    termite07     termite08  2300
148  traj00100    termite08     termite01  1824
149  traj00100    termite08     termite02  1935
150  traj00100    termite08     termite03  1383
151  traj00100    termite08     termite04  1766
152  traj00100    termite08     termite05  1912
153  traj00100    termite08     termite06  1362
154  traj00100    termite08     termite07  2302
155  traj00143    termite01     termite02  1233
156  traj00143    termite01     termite03  1156
157  traj00143    termite02     termite01  1224
158  traj00143    termite02     termite03  2247
159  traj00143    termite03     termite01  1157
160  traj00143    termite03     termite02  2215
161  traj00144    termite01     termite02  1253
162  traj00144    termite01     termite03  1794
163  traj00144    termite01     termite04  1010
164  traj00144    termite02     termite01  1253
165  traj00144    termite02     termite03  8963
166  traj00144    termite02     termite04  4449
167  traj00144    termite03     termite01  1787
168  traj00144    termite03     termite02  8932
169  traj00144    termite03     termite04  8468
170  traj00144    termite04     termite01  1020
171  traj00144    termite04     termite02  4426
172  traj00144    termite04     termite03  8483
173  traj00145    termite01     termite02   656
174  traj00145    termite01     termite03   548
175  traj00145    termite01     termite04   401
176  traj00145    termite01     termite05   308
177  traj00145    termite01     termite06   368
178  traj00145    termite01     termite07   417
179  traj00145    termite01     termite08   566
180  traj00145    termite02     termite01   640
181  traj00145    termite02     termite03   466
182  traj00145    termite02     termite04   555
183  traj00145    termite02     termite05   738
184  traj00145    termite02     termite06   385
185  traj00145    termite02     termite07   439
186  traj00145    termite02     termite08   342
187  traj00145    termite03     termite01   544
188  traj00145    termite03     termite02   472
189  traj00145    termite03     termite04   629
190  traj00145    termite03     termite05   922
191  traj00145    termite03     termite06   594
192  traj00145    termite03     termite07   261
193  traj00145    termite03     termite08   295
194  traj00145    termite04     termite01   395
195  traj00145    termite04     termite02   556
196  traj00145    termite04     termite03   638
197  traj00145    termite04     termite05   618
198  traj00145    termite04     termite06   376
199  traj00145    termite04     termite07   922
200  traj00145    termite04     termite08   549
201  traj00145    termite05     termite01   309
202  traj00145    termite05     termite02   732
203  traj00145    termite05     termite03   926
204  traj00145    termite05     termite04   616
205  traj00145    termite05     termite06   517
206  traj00145    termite05     termite07   581
207  traj00145    termite05     termite08   377
208  traj00145    termite06     termite01   362
209  traj00145    termite06     termite02   386
210  traj00145    termite06     termite03   592
211  traj00145    termite06     termite04   375
212  traj00145    termite06     termite05   515
213  traj00145    termite06     termite07   268
214  traj00145    termite06     termite08   553
215  traj00145    termite07     termite01   423
216  traj00145    termite07     termite02   439
217  traj00145    termite07     termite03   268
218  traj00145    termite07     termite04   923
219  traj00145    termite07     termite05   595
220  traj00145    termite07     termite06   270
221  traj00145    termite07     termite08   396
222  traj00145    termite08     termite01   568
223  traj00145    termite08     termite02   336
224  traj00145    termite08     termite03   298
225  traj00145    termite08     termite04   549
226  traj00145    termite08     termite05   371
227  traj00145    termite08     termite06   553
228  traj00145    termite08     termite07   391
229  traj00013    termite01     termite02    34
230  traj00013    termite01     termite03   254
231  traj00013    termite01     termite04   105
232  traj00013    termite01     termite05    15
233  traj00013    termite01     termite06   420
234  traj00013    termite01     termite07   294
235  traj00013    termite01     termite08    23
236  traj00013    termite01     termite09    24
237  traj00013    termite01     termite10    49
238  traj00013    termite01     termite11    29
239  traj00013    termite01     termite12     1
240  traj00013    termite01     termite13     5
241  traj00013    termite01     termite14  2543
242  traj00013    termite01     termite15    44
243  traj00013    termite01     termite16    19
244  traj00013    termite02     termite01    34
245  traj00013    termite02     termite03   526
246  traj00013    termite02     termite04   535
247  traj00013    termite02     termite05   344
248  traj00013    termite02     termite06   226
249  traj00013    termite02     termite07   378
250  traj00013    termite02     termite08   235
251  traj00013    termite02     termite09   217
252  traj00013    termite02     termite10   410
253  traj00013    termite02     termite11   435
254  traj00013    termite02     termite12   321
255  traj00013    termite02     termite13   500
256  traj00013    termite02     termite14   189
257  traj00013    termite02     termite15   239
258  traj00013    termite02     termite16   427
259  traj00013    termite03     termite01   254
260  traj00013    termite03     termite02   509
261  traj00013    termite03     termite04   508
262  traj00013    termite03     termite05   387
263  traj00013    termite03     termite06   289
264  traj00013    termite03     termite07   193
265  traj00013    termite03     termite08   240
266  traj00013    termite03     termite09   277
267  traj00013    termite03     termite10   262
268  traj00013    termite03     termite11   365
269  traj00013    termite03     termite12   316
270  traj00013    termite03     termite13   387
271  traj00013    termite03     termite14   124
272  traj00013    termite03     termite15   278
273  traj00013    termite03     termite16   522
274  traj00013    termite04     termite01   104
275  traj00013    termite04     termite02   524
276  traj00013    termite04     termite03   494
277  traj00013    termite04     termite05   407
278  traj00013    termite04     termite06   128
279  traj00013    termite04     termite07   497
280  traj00013    termite04     termite08   281
281  traj00013    termite04     termite09   339
282  traj00013    termite04     termite10   312
283  traj00013    termite04     termite11   456
284  traj00013    termite04     termite12   394
285  traj00013    termite04     termite13   257
286  traj00013    termite04     termite14   384
287  traj00013    termite04     termite15   386
288  traj00013    termite04     termite16   278
289  traj00013    termite05     termite01    15
290  traj00013    termite05     termite02   334
291  traj00013    termite05     termite03   363
292  traj00013    termite05     termite04   410
293  traj00013    termite05     termite06   466
294  traj00013    termite05     termite07   436
295  traj00013    termite05     termite08   508
296  traj00013    termite05     termite09   234
297  traj00013    termite05     termite10   383
298  traj00013    termite05     termite11   266
299  traj00013    termite05     termite12   280
300  traj00013    termite05     termite13   398
301  traj00013    termite05     termite14    58
302  traj00013    termite05     termite15   320
303  traj00013    termite05     termite16   358
304  traj00013    termite06     termite01   421
305  traj00013    termite06     termite02   215
306  traj00013    termite06     termite03   275
307  traj00013    termite06     termite04   125
308  traj00013    termite06     termite05   483
309  traj00013    termite06     termite07   303
310  traj00013    termite06     termite08   234
311  traj00013    termite06     termite09   161
312  traj00013    termite06     termite10   183
313  traj00013    termite06     termite11   381
314  traj00013    termite06     termite12   328
315  traj00013    termite06     termite13   345
316  traj00013    termite06     termite14   227
317  traj00013    termite06     termite15   211
318  traj00013    termite06     termite16   259
319  traj00013    termite07     termite01   295
320  traj00013    termite07     termite02   392
321  traj00013    termite07     termite03   183
322  traj00013    termite07     termite04   506
323  traj00013    termite07     termite05   439
324  traj00013    termite07     termite06   293
325  traj00013    termite07     termite08   361
326  traj00013    termite07     termite09   444
327  traj00013    termite07     termite10   242
328  traj00013    termite07     termite11   184
329  traj00013    termite07     termite12   321
330  traj00013    termite07     termite13   421
331  traj00013    termite07     termite14   135
332  traj00013    termite07     termite15   351
333  traj00013    termite07     termite16   364
334  traj00013    termite08     termite01    23
335  traj00013    termite08     termite02   243
336  traj00013    termite08     termite03   228
337  traj00013    termite08     termite04   283
338  traj00013    termite08     termite05   494
339  traj00013    termite08     termite06   210
340  traj00013    termite08     termite07   331
341  traj00013    termite08     termite09   120
342  traj00013    termite08     termite10   395
343  traj00013    termite08     termite11   334
344  traj00013    termite08     termite12   209
345  traj00013    termite08     termite13   543
346  traj00013    termite08     termite14   299
347  traj00013    termite08     termite15   274
348  traj00013    termite08     termite16   293
349  traj00013    termite09     termite01    24
350  traj00013    termite09     termite02   211
351  traj00013    termite09     termite03   278
352  traj00013    termite09     termite04   344
353  traj00013    termite09     termite05   211
354  traj00013    termite09     termite06   171
355  traj00013    termite09     termite07   435
356  traj00013    termite09     termite08   130
357  traj00013    termite09     termite10   244
358  traj00013    termite09     termite11   427
359  traj00013    termite09     termite12   292
360  traj00013    termite09     termite13   410
361  traj00013    termite09     termite14   287
362  traj00013    termite09     termite15   262
363  traj00013    termite09     termite16   200
364  traj00013    termite10     termite01    49
365  traj00013    termite10     termite02   359
366  traj00013    termite10     termite03   251
367  traj00013    termite10     termite04   290
368  traj00013    termite10     termite05   351
369  traj00013    termite10     termite06   189
370  traj00013    termite10     termite07   239
371  traj00013    termite10     termite08   383
372  traj00013    termite10     termite09   222
373  traj00013    termite10     termite11   506
374  traj00013    termite10     termite12   300
375  traj00013    termite10     termite13   376
376  traj00013    termite10     termite14   119
377  traj00013    termite10     termite15   231
378  traj00013    termite10     termite16   265
379  traj00013    termite11     termite01    29
380  traj00013    termite11     termite02   401
381  traj00013    termite11     termite03   361
382  traj00013    termite11     termite04   470
383  traj00013    termite11     termite05   252
384  traj00013    termite11     termite06   400
385  traj00013    termite11     termite07   183
386  traj00013    termite11     termite08   379
387  traj00013    termite11     termite09   442
388  traj00013    termite11     termite10   550
389  traj00013    termite11     termite12   395
390  traj00013    termite11     termite13   183
391  traj00013    termite11     termite14   118
392  traj00013    termite11     termite15   258
393  traj00013    termite11     termite16   357
394  traj00013    termite12     termite01     1
395  traj00013    termite12     termite02   290
396  traj00013    termite12     termite03   319
397  traj00013    termite12     termite04   379
398  traj00013    termite12     termite05   299
399  traj00013    termite12     termite06   314
400  traj00013    termite12     termite07   321
401  traj00013    termite12     termite08   202
402  traj00013    termite12     termite09   271
403  traj00013    termite12     termite10   313
404  traj00013    termite12     termite11   386
405  traj00013    termite12     termite13   304
406  traj00013    termite12     termite14    99
407  traj00013    termite12     termite15   204
408  traj00013    termite12     termite16   517
409  traj00013    termite13     termite01     5
410  traj00013    termite13     termite02   512
411  traj00013    termite13     termite03   405
412  traj00013    termite13     termite04   261
413  traj00013    termite13     termite05   392
414  traj00013    termite13     termite06   374
415  traj00013    termite13     termite07   417
416  traj00013    termite13     termite08   566
417  traj00013    termite13     termite09   404
418  traj00013    termite13     termite10   387
419  traj00013    termite13     termite11   181
420  traj00013    termite13     termite12   286
421  traj00013    termite13     termite14    81
422  traj00013    termite13     termite15   191
423  traj00013    termite13     termite16   214
424  traj00013    termite14     termite01  2543
425  traj00013    termite14     termite02   190
426  traj00013    termite14     termite03   124
427  traj00013    termite14     termite04   391
428  traj00013    termite14     termite05    55
429  traj00013    termite14     termite06   237
430  traj00013    termite14     termite07   137
431  traj00013    termite14     termite08   300
432  traj00013    termite14     termite09   288
433  traj00013    termite14     termite10   120
434  traj00013    termite14     termite11   125
435  traj00013    termite14     termite12    98
436  traj00013    termite14     termite13    81
437  traj00013    termite14     termite15    62
438  traj00013    termite14     termite16   173
439  traj00013    termite15     termite01    44
440  traj00013    termite15     termite02   262
441  traj00013    termite15     termite03   264
442  traj00013    termite15     termite04   376
443  traj00013    termite15     termite05   341
444  traj00013    termite15     termite06   234
445  traj00013    termite15     termite07   354
446  traj00013    termite15     termite08   287
447  traj00013    termite15     termite09   240
448  traj00013    termite15     termite10   247
449  traj00013    termite15     termite11   246
450  traj00013    termite15     termite12   211
451  traj00013    termite15     termite13   210
452  traj00013    termite15     termite14    63
453  traj00013    termite15     termite16   294
454  traj00013    termite16     termite01    18
455  traj00013    termite16     termite02   474
456  traj00013    termite16     termite03   533
457  traj00013    termite16     termite04   280
458  traj00013    termite16     termite05   338
459  traj00013    termite16     termite06   255
460  traj00013    termite16     termite07   349
461  traj00013    termite16     termite08   290
462  traj00013    termite16     termite09   190
463  traj00013    termite16     termite10   269
464  traj00013    termite16     termite11   337
465  traj00013    termite16     termite12   516
466  traj00013    termite16     termite13   224
467  traj00013    termite16     termite14   175
468  traj00013    termite16     termite15   299
469  traj00130    termite01     termite02   951
470  traj00130    termite01     termite03   906
471  traj00130    termite01     termite04   797
472  traj00130    termite01     termite05   668
473  traj00130    termite01     termite06   906
474  traj00130    termite01     termite07   897
475  traj00130    termite01     termite08   736
476  traj00130    termite01     termite09   899
477  traj00130    termite01     termite10  1182
478  traj00130    termite01     termite11   928
479  traj00130    termite01     termite12  1023
480  traj00130    termite01     termite13  1117
481  traj00130    termite01     termite14  1189
482  traj00130    termite01     termite15  1353
483  traj00130    termite01     termite16   876
484  traj00130    termite02     termite01   927
485  traj00130    termite02     termite03   631
486  traj00130    termite02     termite04   646
487  traj00130    termite02     termite05   673
488  traj00130    termite02     termite06   667
489  traj00130    termite02     termite07   925
490  traj00130    termite02     termite08  1107
491  traj00130    termite02     termite09  1040
492  traj00130    termite02     termite10   648
493  traj00130    termite02     termite11   850
494  traj00130    termite02     termite12  1008
495  traj00130    termite02     termite13  1015
496  traj00130    termite02     termite14   630
497  traj00130    termite02     termite15   687
498  traj00130    termite02     termite16  1105
499  traj00130    termite03     termite01   913
500  traj00130    termite03     termite02   659
501  traj00130    termite03     termite04  1269
502  traj00130    termite03     termite05  1026
503  traj00130    termite03     termite06   728
504  traj00130    termite03     termite07   804
505  traj00130    termite03     termite08  1463
506  traj00130    termite03     termite09   793
507  traj00130    termite03     termite10  1081
508  traj00130    termite03     termite11  1007
509  traj00130    termite03     termite12   826
510  traj00130    termite03     termite13   845
511  traj00130    termite03     termite14  1148
512  traj00130    termite03     termite15   880
513  traj00130    termite03     termite16  1090
514  traj00130    termite04     termite01   796
515  traj00130    termite04     termite02   645
516  traj00130    termite04     termite03  1268
517  traj00130    termite04     termite05  1411
518  traj00130    termite04     termite06  1056
519  traj00130    termite04     termite07   810
520  traj00130    termite04     termite08  1694
521  traj00130    termite04     termite09  1685
522  traj00130    termite04     termite10   953
523  traj00130    termite04     termite11  1174
524  traj00130    termite04     termite12   387
525  traj00130    termite04     termite13  1809
526  traj00130    termite04     termite14  1235
527  traj00130    termite04     termite15   645
528  traj00130    termite04     termite16   784
529  traj00130    termite05     termite01   659
530  traj00130    termite05     termite02   683
531  traj00130    termite05     termite03   968
532  traj00130    termite05     termite04  1409
533  traj00130    termite05     termite06  1299
534  traj00130    termite05     termite07   521
535  traj00130    termite05     termite08   736
536  traj00130    termite05     termite09   810
537  traj00130    termite05     termite10   803
538  traj00130    termite05     termite11   801
539  traj00130    termite05     termite12   844
540  traj00130    termite05     termite13   978
541  traj00130    termite05     termite14   945
542  traj00130    termite05     termite15   965
543  traj00130    termite05     termite16   833
544  traj00130    termite06     termite01   879
545  traj00130    termite06     termite02   670
546  traj00130    termite06     termite03   722
547  traj00130    termite06     termite04  1056
548  traj00130    termite06     termite05  1270
549  traj00130    termite06     termite07   824
550  traj00130    termite06     termite08   679
551  traj00130    termite06     termite09   637
552  traj00130    termite06     termite10  1103
553  traj00130    termite06     termite11  1100
554  traj00130    termite06     termite12  1242
555  traj00130    termite06     termite13   887
556  traj00130    termite06     termite14   657
557  traj00130    termite06     termite15   979
558  traj00130    termite06     termite16   601
559  traj00130    termite07     termite01   912
560  traj00130    termite07     termite02   952
561  traj00130    termite07     termite03   784
562  traj00130    termite07     termite04   809
563  traj00130    termite07     termite05   521
564  traj00130    termite07     termite06   822
565  traj00130    termite07     termite08   997
566  traj00130    termite07     termite09   718
567  traj00130    termite07     termite10  1134
568  traj00130    termite07     termite11   933
569  traj00130    termite07     termite12   840
570  traj00130    termite07     termite13   844
571  traj00130    termite07     termite14   835
572  traj00130    termite07     termite15   941
573  traj00130    termite07     termite16   894
574  traj00130    termite08     termite01   727
575  traj00130    termite08     termite02  1092
576  traj00130    termite08     termite03  1415
577  traj00130    termite08     termite04  1692
578  traj00130    termite08     termite05   730
579  traj00130    termite08     termite06   673
580  traj00130    termite08     termite07   998
581  traj00130    termite08     termite09  1017
582  traj00130    termite08     termite10  1081
583  traj00130    termite08     termite11  1082
584  traj00130    termite08     termite12  1129
585  traj00130    termite08     termite13   802
586  traj00130    termite08     termite14  1289
587  traj00130    termite08     termite15  1025
588  traj00130    termite08     termite16   776
589  traj00130    termite09     termite01   894
590  traj00130    termite09     termite02  1045
591  traj00130    termite09     termite03   759
592  traj00130    termite09     termite04  1683
593  traj00130    termite09     termite05   802
594  traj00130    termite09     termite06   654
595  traj00130    termite09     termite07   697
596  traj00130    termite09     termite08  1014
597  traj00130    termite09     termite10  1220
598  traj00130    termite09     termite11   857
599  traj00130    termite09     termite12   888
600  traj00130    termite09     termite13   995
601  traj00130    termite09     termite14  1145
602  traj00130    termite09     termite15   905
603  traj00130    termite09     termite16   752
604  traj00130    termite10     termite01  1163
605  traj00130    termite10     termite02   648
606  traj00130    termite10     termite03  1058
607  traj00130    termite10     termite04   952
608  traj00130    termite10     termite05   768
609  traj00130    termite10     termite06  1093
610  traj00130    termite10     termite07  1137
611  traj00130    termite10     termite08  1108
612  traj00130    termite10     termite09  1217
613  traj00130    termite10     termite11  1011
614  traj00130    termite10     termite12  1053
615  traj00130    termite10     termite13  1325
616  traj00130    termite10     termite14   999
617  traj00130    termite10     termite15  1313
618  traj00130    termite10     termite16  1610
619  traj00130    termite11     termite01   912
620  traj00130    termite11     termite02   864
621  traj00130    termite11     termite03  1018
622  traj00130    termite11     termite04  1181
623  traj00130    termite11     termite05   799
624  traj00130    termite11     termite06  1111
625  traj00130    termite11     termite07   929
626  traj00130    termite11     termite08  1125
627  traj00130    termite11     termite09   862
628  traj00130    termite11     termite10  1030
629  traj00130    termite11     termite12  1025
630  traj00130    termite11     termite13  1094
631  traj00130    termite11     termite14  1367
632  traj00130    termite11     termite15   810
633  traj00130    termite11     termite16  1099
634  traj00130    termite12     termite01  1002
635  traj00130    termite12     termite02   987
636  traj00130    termite12     termite03   817
637  traj00130    termite12     termite04   388
638  traj00130    termite12     termite05   823
639  traj00130    termite12     termite06  1247
640  traj00130    termite12     termite07   839
641  traj00130    termite12     termite08  1127
642  traj00130    termite12     termite09   894
643  traj00130    termite12     termite10  1045
644  traj00130    termite12     termite11  1016
645  traj00130    termite12     termite13   960
646  traj00130    termite12     termite14   753
647  traj00130    termite12     termite15  1139
648  traj00130    termite12     termite16   894
649  traj00130    termite13     termite01  1124
650  traj00130    termite13     termite02  1040
651  traj00130    termite13     termite03   837
652  traj00130    termite13     termite04  1809
653  traj00130    termite13     termite05   980
654  traj00130    termite13     termite06   894
655  traj00130    termite13     termite07   851
656  traj00130    termite13     termite08   825
657  traj00130    termite13     termite09   999
658  traj00130    termite13     termite10  1330
659  traj00130    termite13     termite11  1082
660  traj00130    termite13     termite12   962
661  traj00130    termite13     termite14  1018
662  traj00130    termite13     termite15   872
663  traj00130    termite13     termite16   757
664  traj00130    termite14     termite01  1170
665  traj00130    termite14     termite02   603
666  traj00130    termite14     termite03  1093
667  traj00130    termite14     termite04  1232
668  traj00130    termite14     termite05   912
669  traj00130    termite14     termite06   650
670  traj00130    termite14     termite07   814
671  traj00130    termite14     termite08  1288
672  traj00130    termite14     termite09  1096
673  traj00130    termite14     termite10   956
674  traj00130    termite14     termite11  1343
675  traj00130    termite14     termite12   755
676  traj00130    termite14     termite13   986
677  traj00130    termite14     termite15  1177
678  traj00130    termite14     termite16   862
679  traj00130    termite15     termite01  1363
680  traj00130    termite15     termite02   682
681  traj00130    termite15     termite03   882
682  traj00130    termite15     termite04   646
683  traj00130    termite15     termite05   966
684  traj00130    termite15     termite06   975
685  traj00130    termite15     termite07   920
686  traj00130    termite15     termite08  1038
687  traj00130    termite15     termite09   879
688  traj00130    termite15     termite10  1329
689  traj00130    termite15     termite11   775
690  traj00130    termite15     termite12  1142
691  traj00130    termite15     termite13   880
692  traj00130    termite15     termite14  1196
693  traj00130    termite15     termite16  1075
694  traj00130    termite16     termite01   861
695  traj00130    termite16     termite02  1108
696  traj00130    termite16     termite03  1054
697  traj00130    termite16     termite04   784
698  traj00130    termite16     termite05   804
699  traj00130    termite16     termite06   612
700  traj00130    termite16     termite07   883
701  traj00130    termite16     termite08   783
702  traj00130    termite16     termite09   761
703  traj00130    termite16     termite10  1583
704  traj00130    termite16     termite11  1059
705  traj00130    termite16     termite12   906
706  traj00130    termite16     termite13   743
707  traj00130    termite16     termite14   875
708  traj00130    termite16     termite15  1071
\end{Soutput}
\end{Schunk}

\end{document}
